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^ Abstract 

D . We consider the effect of warping on the distribution of type IIB flux vacua constructed 



with Calabi-Yau orientifolds. We derive an analytical form of the distribution that incor- 

CN ■ porates warping and find close agreement with the results of a Monte Carlo enumeration of 

(-\| I vacua. Compared with calculations that neglect warping, we find that for any finite volume 

I compactification, the density of vacua is highly diluted in close proximity to the conifold 

. point, with a steep drop-off within a critical distance. 
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1 Introduction 

^ I Complex structure moduli in type IIB string theory are stabilized by turning on fluxes, and in 
H ! 

. certain parts of the theory's moduli space the fluxes lead to large warping effects. These effects are 
essential for a detailed understanding of dynamics in the string theory landscape [1]. Tunneling 
between flux vacua involves the nucleation of a brane carrying appropriate charges, but such 
events appear to be favored in conflgurations of the Calabi-Yau geometry where particular cycles 
are small, and hence, warping due to fluxes through such cycles is large. 

Aside from dynamics, the distribution of vacua for such models is of interest. The Bousso- 
Polchinski model of the string landscape [2j suggests that with a sufficient number of fluxes, 
one should expect vacuum energies that are sufficiently flnely spaced to ensure that some vacua 
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have cosmological constants in rough agreement with our own. In [31 II], the authors developed 
a convenient framework for carrying out such analyses in the context of type IIB string theory 
compactifications. The theoretical vacuum distributions for certain simple Calabi-Yau compacti- 
fications have also been supported by numerical studies [5]. 

A quite general result of these studies is that vacua appear to accumulate around the conifold 
point in the complex structure moduli space. In fact, the density of these vacua diverges logarith- 
mically. However, in light of the fact that warping becomes strong precisely near the conifold point 
for any finite volume Calabi-Yau compactification, the natural question arises of what effect — if 
any — warping may have on the distribution of vacua. 

A related point is whether the vacuum density is well-captured by the simpler to compute 
index density. In the absence of warping, the index is simpler to compute since it is similar to 
the Chern class of the moduli space, whereas there is no straightforward geometric quantity that 
corresponds to the vacuum count. As we shall see, when warping is included the overall agreement 
between number and index densities continues to hold, but the topological nature of the latter 
becomes more complicated. 

In section [3] we review the general framework for deriving theoretical distributions of vacua 
for unwarped Calabi-Yau compactifications as originally laid out in [4J. We then explain how to 
modify this construction to derive the warped version of the number and index densities. The 
methods are then used to explicitly compute the densities in the vicinity of a conifold point. In 
section H] we numerically generate near-conifold distributions of vacua and compare these to the 
theoretical distributions of section [21 

2 Background 

Type IIB string theory compactified on a Calabi-Yau manifold yields scalar fields, known as moduli, 
in the low energy supergravity theory. These moduli are related to geometrical parameters of the 
internal Calabi-Yau manifold and can in certain models number in the hundreds. Unfortunately, 
these moduli appear as massless fields without any potential governing their dynamics, rendering 
the physics unrealistic. Fortunately string theory contains other ingredients with the capacity to 
resolve this problem. In particular, these scalar fields can be stabilized by turning on various p-form 
fiuxes in the internal manifold, a procedure that generates the Gukov-Vafa-Witten superpotential: 



Here, Q3 is the holomorphic (3, 0) form defined on the Calabi-Yau, G3 = F3 — tH^ is the type IIB 
3-form field strength, r is the axio-dilaton, and z denotes the set of complex moduli mentioned 
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above upon which the holomorphic three form depends. Given this superpotential, the scalar 
potential for the moduli becomes: 

V{z,r) = e^/^^ (^K'^-^DaWDW-^m'^ (2) 

where the sum runs over the complex moduli = 1,2, ... , n), with n = h^y, as well as the 
axio-dilaton = 0). Here, the covariant derivative acts as = da + where Ka is the 
derivative of the Kahler potential with repsect to the a^^ complex modulus or r. This ensures 
that DW transforms in the same way as W itself under a Kahler transformation so that the 
physical potential V{z,t) is invariant under Kahler transformations. Supersymmetric minima of 
the potential V occur at points in moduli space where DaW = with a running over all of the 
moduli and the axio-dilaton. In general, the potential has many minima, each of which represents 
a stable low energy configuration of the internal Calabi-Yau. These configurations arise from the 
large number of discrete fluxes that can thread through the Calabi-Yau's various 3-cycles. It is 
thus natural to explore this large landscape of flux vacua using statistical methods, as was first 
done in [31 S], as we now briefly review. 



3 Analytical Distributions 
3.1 Counting the vacua 

Here we will review the derivation of the index density given by Douglas and Denef in |3], focusing 
on areas where our analysis, including the effects of warping, will differ. We will restrict attention 
to vacua that satisfy DaW = for all complex moduli and the axio-dilaton. The strategy is to 
consider these equations as constraints on the choice of fluxes and otherwise, simply allow the 
fluxes to scan. First, assume that fluxes are fixed and consider the function on moduli space given 

b>0 

6^^+^{DWiz)) = 6iDoWiz)) . . . 5iDnWiz))6iDoWiz)) . . . 5(D„W^(^)). (3) 

Clearly this provides support only at the locations of the vacua. However, as written each vacuum 
does not contribute with the same weight. To see this, rewrite equation ([3]) as a sum of delta 
functions which explicitly spike at the locations of the minima: 

= (4) 



^Our conventions for the delta functions and integration measures depending on a complex variable z are given 
by S'^{z) — S(Re z)5{lm z), and (Pz — d(Re z)d{lm z). 
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Here the determinant arises from expanding the delta functions near each minimum in much the 
same way as 6{f{x)) = ^ 6{x — X2eTo)/\f'{x)\, and is of the {2n + 2) x [2n + 2) matrix 



daD.W daDkW\ 

where we let a, h range over the n moduli as well as the axio-dilaton. Note that the partial 
deriviatives in the matrix above can be replaced by covariant derivatives at the vacua since there 
the conditions DaW = render the two expressions equivalent. If we then integrate this over the 
moduli space we find contributions from each vacuum associated with a fixed set of fluxes with 
weight I det D'^W\~^. Since this value is not constant over the moduli space, the result will not 
reflect the number of vacua. To count the vacua, we must compensate by integrating over the 
delta-functions appropriately weighted: 

J d^''zd^T6^''+\DW{z))\detD^W\. (6) 

This expression defines the vacuum count for a given set of fluxes. Another useful quantity 
considered in [1] is the index, which involves dropping the absolute values around the determinant 
of the fermion mass matrix: 

J d^'^zd^T d^'^+^DWiz)) det D^W. (7) 

This integral then counts the number of positive vacua minus the number of negative vacua, where 
parity is given by the sign of the determinant of the matrix in equation ([5]). To count all vacua, 
we must also sum over fluxes subject to the tadpole cancellation condition 

F3AH,<L,. (8) 

CY 

Here is the maximum possible value for L. It will turn out to be useful to lift this discussion 
to F-theory where we consider our manifold M as an elliptically flbered Calabi-Yau 4-fold, whose 
base consists of the original 3-fold and flbers are given by the auxiliary 2-torus whose period is 
given by the axio-dilaton r. We decompose the holomorphic 4-form: 

^4 = 1^1 A fig, (9) 

where Qi is the holomorphic one form on the two torus parameterizing the axio-dilaton, and 
is the usual holomorphic three form on the Calabi-Yau. In particular, if we consider the two 
one-cylces A and B on the torus, we can deflne the two one forms a and /3 dual to the cycles A 
and B such that /_4 7 = « A 7 and /g7 = Jrp2 /3 A 7 for all closed one forms 7. Then, as long 
as we deflne our holomorphic one-form Qi as 

ni=a-T^, (10) 
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we will have r = J^^i/ /g^i as we want for the complex structure of the torus. Furthermore, if 
we define a flux four form as 6*4 = /3 A F3 — a A ifs, we can write the tadpole condition as 

1 



G^AGi = - a A (3 F3AH3 (11) 

M Jc'Y'j, 

If we normalize the F-theory torus volume so that « A /3 = —1, this exactly reproduces 
the tadpole condition in the type JIB picture. With K = dim H^y^ we've lumped the 2K fluxes 
Fq, . . . , Fk-i, Hq, . . . , Hk-i into the 2K components of G4. Also note that with this definition of 
the flux four form, we can write the usual type JIB superpotential as 



W= n^AG^. (12) 

J M 

Let's choose a particular basis of three forms on the CY3 {Sj}, and denote the intersection 
form in this basis as Qij so that 

^ S,AS, = Qi,- (13) 



We can extend this basis to M. by wedging it with the one forms a and (3. In this basis, we denote 
the components of the field strength G4, by A^^ with a = 0, 1, . . . 2K — 1, and the intersection form 
in the full 4 (complex) dimensional space by rjab- Then, the tadpole condition in equation ([8]) can 
be written in terms of the components of the two fluxes {F = F^Y^i and H = iJ*Sj) as 

L = ^N'^VabN' = F'Q,,H^ < (14) 

We should then sum only over the fluxes that satisfy this inequality. In particular we can imagine 
summing over all fluxes while including a step function. 



Index = ^e{L,-L) j d^'^zd^r 5^''+\DW{z)) det D^W 



(15) 



Fluxes 

We can write the step function as an integral over a delta functio: 



e{L,-L) = / 6{L-L)dL (16) 
J —00 

yielding 

Index = ^ f dl f d'^''zd\ S{L- 1)6^''+'^ {DW{z)) det D^W (17) 

Fluxes 

By treating the fluxes Nq, . . . , N2K-1 as continuously varying parameters, we can approximate this 
sum by an integral. 

Index = [ * dl [ d^^N [ d^^'zd^r 5{L - 1)5^''-^^ {DW{z)) det D^W (18) 



^Note that in [4j the step-function is expressed in terms of a contour integral over an exponential e"^*. Our 
expression in terms of a delta function proves to be more useful for the analysis incorporating warping effects. 
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It is natural to define the index density in moduli (and axio-dilaton) space by 

^lI{z,T)= I ' dl I d^^N 5{L~L)5'''^+''{DW{z))deiD''W (19) 

J — oo J 

Upon integrating over r, this will then equal the total index. We now rewrite this index 
density in terms of geometric properties of the moduli space. A first step in doing this is to change 
basis from {a A S^, /3 A S^} to the set of linearly independent four forms {f24, Da^^-, D^DiVti} U 
{c.c} where a ranges over the complex moduli as well as the axio-dilaton while i ranges only over 
the moduli. This proposed basis consists of 4(n + 1) elements where n denotes the number of 
complex moduli in our theory, which agrees with the 2K elements of the original basis. This new 
basis satisfies 



= e-^("'") (20) 



M 



[ Dan,AD-,n, = -e-^(^'^)ir,5 (21) 
Jm 

[ DoDiQ^ A D-oD-pi = e-^^^^''^ K^fKfj, (22) 
Jm 

with all other combinations vanishing. By rescaling all of our basis elements by the factor e^^'^'^^^'^, 
the new basis won't have any of the extra exponentials in their inner products: 

[ e^(^'^)/2^4 A e^(^'^)/2^4 = 1 (23) 
Jm 

f e^("'^)/2^ifi4 A e^("'")/2Djn4 = -Kfj (24) 
Jm 

[ e^("'")/2^oAf^4 A e^("'")/2Dor)j04 = K^fKfj, (25) 
Jm 

And becuase of the properties of the covariant derivative, we can accomplish these changes by 
rescaling the holomorphic 4-form by this same factor: — )■ e^'^'^'^^^'^ri^. For notational simplicity 
we will redefine ^4 to represent this rescaled versio When we want to explicitly refer to the 
actual holomorphic 4-form, we will denote it as ^4: 

= e^(^'^)/2^4 (26) 



The covariant derivative Da = da + Ka, is the Hermitiaii metric connection acting on sections of the complex 
hne bundle i, where is a section oi H ® with H the Hodge bundle and L is a line bundle whose first Chern 
class is the Kahler form on the 3-fold's moduli space. The expression J A Ha provides a metric on L from 
which the metric connection then follows. When acting on sections of other, related bundles, the Hermitian metric 
connection must be appropriately modified. 
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Finally, we can consider the set B = {^4^, D^^^, DqDjQ4} U {c.c.} where = e^Da, and 
the vielbeins satisfy e^e^^K^i = 6^^, as usual. The notation is consistent assuming a suitably 
defined spin-connection (see appendix IA.1|) . Our new basis is orthonormal: 



n^An^ = 1 (27) 

M 

DAVtiADsVLi = -Sab (28) 

M 

[ DoDin^AD-oDj-n^ = Sjj, (29) 
Jm 

The 4-form flux G4 in the new basis is given by 

= Xn^ - Y^Da^4 + Z^DoDin^ + c.c. (30) 

with X, Y^, Z^, X, y"^, Z^ being the coefficients of G4 in this basis. Note that G4 does not depend 
on the complex structure or axio-dilaton, which implies that the coefficients X, Y^, Z^ , . . . depend 
on and r in a way that precisely cancels the dependences arising from and its derivatives. 
Since G4 doesn't depend on the complex structure of the Calabi-Yau or the axio-dilaton, we can 
relate these coefficients to various combinations of derivatives acting on the superpotential. In 
particular 



w = 


J 1^4 A ^4 = X 


(31) 


DaW = 


J Da^a AGi = YA 


(32) 


DoDqW = 





(33) 


DoDiW = 


j DoDjn^ AGi = Zi 


(34) 


DiDjW = 


J DjDjQ^ A G4 = J'ijkZ'' 


(35) 


DjDjW = 




(36) 


DqDoW = 


X 


(37) 


DqDiW = 





(38) 



where the computations establishing these relations are provided in the appendix. Note that we 
have defined the coefficients Fijk = i Jcy^^ ^ DiDjDk^^ = i J^y^^ ^ didjdx^s- Also, note 
that W denotes the rescaled superpotential; when we want to explicitly refer to the original one, 
we will once again place a hat on top of it (W). We can then rewrite our expressions in terms of 
these new functions on moduli space, and in particular have for the tadpole condition 

L = -Nr]N =- I GiAGi = \X\^ - \Y\^ + \Z\^, (39) 
2 2 ./ 



where = Y^Y^d^j^, etc. The index density then becomes 

2n+2v^2n7 T I j^. „| r/r _ I TL^|2 , IV|2 _ I 7|2u2n+2/ 



^j{z,t)= I dL I d'Xd'''+'Yd'''ZJ\detg\6{L-\X\' + \Y\'-\Z\')5'''+\YA)\X\ 

X det 




J^ijkZ^ 
Z^ X6ij - ^ 



(40) 



X 

Here J is the Jacobian obtained in changing variables from A^^ to X, Ya, Zj, which we will 
determine explicitly below. We have included an additional factor of | det g\ which comes from 
transforming both the delta functions and the determinant to the new variables, and note that 
factors of cancel between the delta functions and the determinant. 

Let's now compute the Jacobian \J\. In the original basis, the components of G4 were given 
by Na- We can now write the Na in the new basis 

N = r]~^ (XU - Y'^DaU + Z^DqDjU + c.c.) (41) 

Here the Hs are the periods of the rescaled holomorphic four form and are related to the usual 
ones by a factor of e^'^^.We can see from this expression that the change of basis is achieved by 
the application of the matrix M = 77"^ (H, —D^n, Dq-D/H, c.c). If we use the convention that 
d'^z = ^dz A dz, we find that the appropriate Jacobian is 

J = 22("+i)|detM| =4"+^|detr/|-^/2|^etM^r/M|^/2 (42) 

We have M^rjM = diag(l, — 1„, 1, — 1„), which follows from our choice of an orthonor- 
mal basis of 4-forms. This implies that the Jacobian is given by 

J = r+'\ det ri\-'/\ (43) 

The final expression is then (after explicitly integrating over Ya), 

fii{z,T) =r+^\ det r]\-^/^ f * dL I d^ Xd^''Z\ det g\6{L - \X\^ -\Z\^)\X\'^ 

J —00 J 

X det I ^ 7 7- (44) 

We can explicitly integrate over the phases, leaving only integrals over the magnitudes |X| 
and \Z\, showing that the tadpole delta function fixes the region of integration to lie on a circle 
of radius \/Z in the |X|, \ Z\ plane. There is therefore no need to integrate over negative Ls, and 
furthermore the remaining finite integral can be evaluated. Following this approach, one can show 
that the index density has a nice geometrical interpretation [1]: 

Hj{z,t) = det{R + uI), (45) 
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where R is the curvature two form on the moduh space and u is the Kahler form. For the case 
of one complex modulus (and the axio-dilaton), this reduces to yU/ = — 7r^| det r/l^^/^Wo A -Ri where 
Wo is the Kahler form on the axio-dilaton side while Ri is the curvature form on the moduli space 
side. In order to obtain this, one must use a relationship between the Kahler and curvature forms 
on the axio-dilaton moduli space: Rq = —2uo. 



3.2 Incorporating warping 

A full treatment of warped Calabi-Yau geometry involves using the machinery of generalized 
complex geometry [HI [71 18|. However, a rough method that produces the appropriate functional 
behavior induced by warping near the conifold will suffice for our purposes. This behavior can be 
derived by taking the warped Kahler potential to be approximated by [H] 

e-^ = f e-^^nACl^ f nACl+ f (l + ^^^nACl, (46) 



J J Bulk 7 Conifold \ ^ / 

where e~^^ = 1 + e~^^°/c is the warp factor, with e~'^^° capturing the significant warping at the 
conifold while c is a constant related to the overall volume of the Calabi-Yau manifold. In general, 
we will use tildes to denote quantities that include warp corrections. 

The warp-corrected Kahler metric has been shown to have the near-conifold form [H [T21 [131 [S] 



where ^ is the local coordinate around the conifold point, k = lim^^^p e^^^'^^ and Ki is a constant 
(to leading order) associated with the Kahler metric's expansion around the conifold. The hatted 
quantity in the rightmost expression corresponds to the warp correction to the original, unwarped 
Kahler metric. The constant Cw is on the order of the inverse volume of the Calabi-Yavo, capturing 
the suppression of the warping effects at large volume. 

Given the form of the Kahler metric near the conifold (H7|) . we find that up to shifts by functions 
holomorphic and antiholomorphic in ^, we have 

K ^ K + 9C^\^\^/'' = K + K, (48) 
^ K^ + 3C^^ = K^ + K^. (49) 

To take warping into account in computing the vacuum count and index, we follow the basic 
logic of section |3TT] with incorporating various necessary modifications. First, we continue to define 
quantities such as X, Y, and Z without making any reference to the warping. This means that the 



■^In fact, it goes approximately like Vqy^^ : since it is related to the universal Kahler modulus zero mode. 
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logic for converting the step function 9{L^ — L) into an integral is unchanged. What does change 
are the expressions within the delta-functions and the determinant of the fermion mass matrix. 
In particular, the positions of the vacua are now determined by the conditions DaW + KaW — 0, 
where the second term is the correction due to warping. Thus, the delta-functions must now read 

5^^+\Ya + KaX), 

and the quantities appearing in the fermion mass matrix now have to incorporate warp corrections: 
i^A + Ka^ (^Db + W . Note that at a vacuum we have the equivalence dA^DsW + KbW) = 

(^Da + KaJ (^Db + Kb^ W, and in general we will make use of similar equivalences in what 
follows. We have: 

Do{Dj + Kj)W = Zj, 

(di + ki^ [dj + /?,;) w = Tukz"^ + kijX + kjYj, 



{di + Ao) (dj ^k^w = [Sij + kij) X. 



Upon integrating over the Ya we find the index density 




d?X (f"" Z g\5 [L - a\X\^ - \Z\'^] \X\ 



where a — 1 — KjK , Hij = 6jj + Kjj, and ajj — Kjj — KjKj. 

In order to compute this density, it proves helpful to consider the special case of one complex 
modulus as well as the axio-dilaton. In this particular case, we obtain the expression 

At7 0c j"^* dL j d''Xd''''Z\detg\5{L-a\X\^ -\Z\''^ {{Z^" + {ii^ -\a\^)\X\' - {2ii+\T\^)\X\^\Z\^) 

(50) 

Note, that we have eliminated a few terms that will integrate to zero because they depend 
explicitly on the phases of X, Z. Far from the conifold, a approaches 1 since the warping corrections 
can then be neglected. However, when one moves toward the conifold, a gets progressively smaller 
until at some critical value it equals zero, and then the warping correction drives a negative. As 
long as a is positive, the tadpole delta function fixes the range of integration so that |X| and \Z\ 
lie on a finite ellipse. Upon computing the integral, one therefore obtains a finite value for the 
index. However, when a goes to zero, this ellipse becomes increasingly stretched until for a = 0, 
the range of integration for \X\ becomes unconstrained. At this point, the integral above for the 
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index density diverges. Then, as a goes negative, this divergence persists as the elhpse turns into 
a hyperbola. Naively, this suggests an infinite number of vacua within a finite disk surrounding 
the conifold point. However, the more careful analysis taking account of finite fiuxes that we carry 
out below yields a finite result. 



3.3 Finite fiuxes 

One major difference between the analysis above and numerical simulations is the range of fiuxes. 
In numerical simulations fiuxes are necessarily kept within a finite range, while in the derivation 
above, arbitrarily large fiuxes were included. To derive a theoretical distribution that mirrors the 
effects seen in numerical studies, it is best to include a bound on the fiuxes in the analysis. This 
complicates the final expression for the theoretical distribution but, of course, the finite bound 
on fiuxes is physically well motivated since the supergravity approximation breaks down for large 
enough fluxes. In the absence of warping, the finite range of fiuxes does not lead to dramatic 
differences from naively taking the bound to infinity, but as we will see, this limit is more involved 
when warping is includedlf] 

Suppose that we bound our fiuxes by the range Ni G [—A, A]. The Na and X,Y,Z variables 
are related by 

X = Nalla (51) 
Ya = NaDAUa (52) 
Zj = NaDoDjUa (53) 

Here the H's are the periods of the rescaled holomorphic form, as before. We would thus expect 
the ranges on X, Y, Z to be moduli dependent. Let's separate the phase and magnitude of X, Y, Z. 
Although in principle, the ranges of the phases may have a complicated dependence on both the 
moduli and the magnitudes |X|, \ Z\, we will neglect this subtlety and suppose that they range 
over the usual [0, 27r]. As a result, we can easily integrate these variables out, leaving us with the 
integrals over the magnitudes. We would expect to have these range over the values 

(54) 
(55) 
(56) 



|X| 


G 


[0,A/x(e,r)] 


\Ya\ 


G 


[0,A/y(e,r)] 


\Zi\ 


G 


[0,A^(e,r)] 



'^By way of comparison, if we express the tadpole condition in the manner of [4J, it leads to an integral over 
a Gaussian-like exponential factor e^^''^/^ = e^'^' +1^' ^'■^l , a damping term in the absence of warping due to 
the SUSY conditions 6{Ya)- Warping modifies these conditions to Ya + KaX — 0, and so the argument of the 
exponential is not negative definite, requiring an additional regulator bounding the fluxes. 
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for particular fx, /y, and fz- Let's consider fx- The largest value that \X\ will take corresponds 
to the fluxes Ng, taking one of their two extreme values of ±A; which of two possibilities maximizes 
\X\ depends on the near conifold behavior of the periods. We must choose the eight signs for the 
eight fluxes Na in such a way that we maximize the expression 



max 



E±na(e,r) 



(57) 



Since the periods are all flnite in the near conifold limit, the ^ dependence decouples. However, 
the value for fx will still be r dependent. As far as fy and fz are concerned, the idea is the same 
except for the fact that the ^ dependence can't be neglected due to logarithmic divergences. In 
particular, we flnd that 



fx = fxir) 

fy = /^(r)|l-/?(r)log(0| 
fz = /l(r)|l-/|(r)log(0| 



(58) 
(59) 
(60) 



Now consider a flxed point in moduli space ^ as well as a fixed value for r. Then, the upper 
limits on these integrals will involve particular constants multiplying the fiux cutoff A. Integrating 
over the variables Yq in the expression for the index density, the delta function S (Yq) fixes Yq — 0, 
leaving us with 



Hi oc 



X 



dL 

-oo Jo 



\X\d\X\ 



\Y\d\Y\ / \Z\d\Z\\detg\S{L - a\X\'^ - \Z\'^) 
Jo 

(Fi + K^X) {\Z\' + - |a|2) \X\* - (2/. + \:F\') \X\'\Z\') (61) 



The remaining delta function constraints come from the tadpole condition and the supersymmetry 
condition D^W + K^W — 0, equivalent to Yi + K^X — 0. Satisfying these constraints will place 
complicated restrictions on the upper and lower bounds of the remaining integrals. Let's first 
examine the region of integration imposed by the supersymmetry constraint: 

• When |X| is at its lower bound of 0, the constraint is trivial to satisfy. Thus, the lower 
bound of |X| is unchanged. 



K^X 



> /yA. At 



K^X must vanish. We thus 



However, when \X\ > 0, there will be points in the moduli space where 
such points, the delta function imposing the constraint Yi = 
see that the upper bound of |X| is restricted in such cases to /yA/ . Solving the delta 
function constraint for Yi requires that the upper bound of the \X\ integral be taken to be 
|X|A = min (A/x,A/y/|^^|). 



12 



Note that in our scheme for bounding the fluxes, the upper hmits of integration for |X|, \ Y\ and 
\Z\ all scale with the cutoff A in the same way. So, simply taking the limit as A — )■ oo won't 
affect the analysis. From our scheme's perspective, it is only in the strictly infinite case where 
the naive divergence reappears as discussed at the end of section 13.21 (One could imagine more 
complicated schemes for bounding the fluxes, treating X, Y, and Z independently, allowing for a 
set of continuous limits that recover the divergent results of the naive approach. Such a scheme 
would increase the difficulty of relating the numerical and theoretical analyses, as investigated in 
unwarped case in |3l IH [5]. 

Given the new limits of integration on we can freely integrate out the delta function fixing 
the value of Yi. 

_ cI^Ia rfz^ , X 

dL \ \X\d\X\ \Z\d\Z\\detg\6(L-a\X\'^ -\Z\'^] 

X {\Z\'+{fi'-\a\')\X\'-{2f,+ \T\')\X\'\Z\') (62) 

To simplify our notation, let's change variables to u = |Xp and v = \Z\^ . The density can then 
be written as 

fx I (X dL du dv\ detg\ 5 iL — au — vj (v + (/x^ — |cr|^) — (2fi + jJ-'P) uv) (63) 

J-oo Jo Jo ^ 

where ua = \X\l = min ^AVl, A^y/ 

It's useful to consider the two cases a > and a < 0, separately. 
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« (c) « (d) 

Figure 1: Various possible regions of integration for a > 0. In (a) the L^, > La, where La = 
auA + /|A^ and u\ = min (^A^fj., K^fy/K'^ so the region is cut off at ua- In (b) /|A^ < < Lj^. 
In (c), L* < /|A^, and ma < L^/a. Finally, (d) shows a region where L* < /|A^ and ma > L^ja. 

3.3.1 The case a > 

The delta function in (163|) arising from the tadpole condition is 5 ^L — — f ^ This constrains 
the value of f to be L — aw, as well as constraining the region of integration on the L/w-plane. 
Let's first determine the lower bounds on L and u: 

• The variables u, v are positive or possibly zero and since a > 0, then L > 0, fixing the lower 
bound of for the L integral. 
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• Let = /|A^ be the upper bound on v. If v^p < L, then the delta function forces the 
lower bound on u to be (^L — fz-^^^ /a. However, if v^p > L then the lower bound on u is 

0. So in general we let the lower bound on u be u^^^^ = max (0, — ^ — j . 
Now for the upper bounds: 

• If < auA + fz^^ 1 then it remains the upper bound for the L integral. If on the contrary, 
the inequality runs the other way, > oma + fz^^i then the constraint L — au — v cannot 
be satisfied everywhere along the range < L < L,,, truncating this range to < L < 
otUK + instead. So the upper bound on the L integral is 



L+ = min (L*,L 



A 



where La = aup^ + /|A . 

• If, at a fixed L, we had a-UA > L, then since the lower bound on v is 0, this places an upper 
bound on u of L/a. If the inequality is reversed, then the upper bound on u is u\. So, in 
general the upper bound on u is u^^ = min (^u\, L/a^. 

Various possible regions of integration in the L/u-plane are illustrated in figured] 
Using the bounds described above and integrating over v yields 

/it oc / dL du\ detg\ { i L — au] + (3u^ + 'ju { L — au] ] 

where and /3 = /x^ — |(Tp and 7 = — 2/i — Then, expanding everything out and integrating 

over u, we obtain 

/^l « ^ dL\ detgllz^ («+ - m+^J + ^L^L )^ - 



where the L dependence of u^p and Udown has been suppressed in the last line. 

In order to integrate over L, we must separate the integral above into two parts since u^^ and 
""down different functions of L. Let be the portion of the integral involving terms containing 
powers of and J^^^^ be the portion of the integral containing Mdown- ^ot^ that we remove the 
|det (/I factor from these integrals. Focusing first on we see that foiOO < L/a < u^ (L*), we 



^Note that we could have considered it+p (i^p) mstead of it+p (i*) as the upper part of the interval. However, 
recall that L^^ is the smaller of either L* or Lj^. If > La, u+p (i^p) = u^p {L\) = u\ since from the definition 
of _La, the intequality u\ < L\/a always holds. 
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can replace instances of (^Lj in the integral with L/a, while if u^^ (L*) < L/a, then u^^ — u\, 
which is independent of L. So J^^+ splits into integrals over the two regions: 



rLtp 
J aut„(Lt 



"Lip / ^, r)^, ^,2 

I 



+ / dLlLW + ^^^Lul + ^^^^-^ul) (64) 



Integrating yields 



up 



a 2a2 3a3 y 4 



+ ^ ^ ^ + ((^up) - « Kp (L,)) ) ul (65) 

+ ^±1^ (L+ - a< (L.)) ul (66) 

For the integral J^^^^, we consider the regions < L < /|A and /|A < L. In the first 
case, M^o^n = in which case this entire portion of the integral vanishes, while in the second, 
^down ^ (-^ ~ fz^) 1^- ^* < /i^^ ^^^^ entirety of J^+^^ = 0, so we have 



2q;2 3o;3 



Integrating yields 



^down - ^ /zA ) (^/zA + ^^^^ J 

- /^A ^-^ + /^A (^^-^ + — j (L„J 



3'^ \a q;2 q;3 y/ ^ V 4a Sa^ 12^3 

16 



Notice that when one ignores warping and the finite fiuxes, a = 1,K^ = 0, and A — t- oo, 
implying /3 = 1, 7 = —2 — l^f, u^^ {L^) = L^,, and = L*. In this case, we must go back to 
the expression and note that the second integral in that expression vanishes since the lower 
and upper bound of integration are both L*. Furthermore the integral ^^^^ vanishes due to the 
^-function prefactor. The index density in the unwarped case is thus 

f^i U,^) = Idet^l ^+ = Idet^l — ( 1 = |det c/| — (2-|J^| ) (67) 

This precise combination gives us the curvature tensor as argued in [3l H]. So, our expression 
reduces to the correct form in the unwarped, infinite fiux case. We now turn our attention to the 
case where a < 0. 



3.3.2 The case a< 

Once again, we first establish the lower and upper bounds on L and u: 

• Suppose we are at the lower bound on v, namely v = 0. In this case, the tadpole constraint 
L — au — f = tells us that the lower bound attained by L is L'^own ~ Q^'^A- Note that this 
is negative. 

• Consider some fixed L < fz-^^', If -^v > 0, then there is always a v = L to cancel it, and the 
lower bound for u in this case is 0. However, if L < 0, the fact that f > implies that for 
the constraint to hold, we need the lower bound for u to he L/a. So in general, the lower 
bound for u is tiJown ~ (o, L/a^. 

• By similar reasoning to the previous case, if < fz-^"^, then it may remain the upper bound 
on L. However, if L^, > fz^^i then the upper bound on L becomes fz^'^- So in general, the 
upper bound on L is = min(L^,, /|A^). 

• Consider again a fixed L, and suppose v is at its upper bound of fz-^"^- If C(U\ > L — fz-^'^, 
then the upper bound of u must be truncated to (^L — /|A^^ /a. Otherwise, if au\ < L — 

fz-^^, then the upper bound on u remains ua- In general then, u'^ = min (^ua, (^L — /|A^ j / a 

Given these bounds on u and L, we may now integrate over v, eliminating the tadpole delta 
function to get 

Hj (X dL _ c/m I det^f] (L^ + (7 - 2a)LM + (a^ + /3 - 07) (68) 

Jl- Ju-(l) ^ ' 



down 
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Carrying out the u integration yields 



f^i 0^ dL\ detg\ (l^ {u^^ - u^^^^) + L f^-^l ((«up)' - Kown)') 

down \ 



where we have suppressed the L dependence of u~p, and u^^^^^. 

As before, spht the integral into two parts, J^^" and J^^own' involving just the u~p and u^^wn 
parts, respectively. To compute we consider two cases: 

• Suppose < La, where we recall La — auA + fz^^- Note that since a < 0, we have that 
Lk < fz^^: thus, L~p = in this case. We also see that — u^, and so in this case, 
the integral is simply 



up 



,7/72 , J [ l-'2a \ 2 , g' + Z^-aT 3 
dL \ L ua + L \ — - — J Ua H u\ 

"^down 



• Suppose that > La- In this case, for L < L\, u^^ — u\ as before, but when L > L\ we 
have u~p — — /|A^^ /a. So the integral splits into two parts 



3 



"^down \ 

La \ 



2 



a I \ z / \ a 

3 ■ 



a 



These two expressions can be joined if we introduce L^^id — min {L^, La): 

/Lmid ~/'y — 2a\ 

—'down ^ ^ ^ 



L'^/A + L { ^ — ) ul + ^ '-ul 



a I \ z / \ a 

3 ■ 



^ a^ + P-a^ i L-flK' ^^^^ 



a 
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The integral in the second hne above vanishes if L^i^, — L*, since in that case L^^ also is L*. 
Carrying out the integral yields (after plugging in L'^^^^ — au\) 

^up = Y ( ^mid -a Uj^j+ ^ Uj^ [ ^mid -aU^] + 5 ^^A (i-mid " OIUa) 



The integral J^^own vanishes when L > since in that case M^own ^ 0- Thus, the only region that 
contributes is where L'^own < L < 0, in which u^^^^ = L/a. We have, 



-f 

Jl- 



^down - / dL[-L + -^-^L + — L 



which gives 



1 /I 7 - 2a + /3 - a7\ 



where we have again used L^^^^ — au\. 
The full index density is thus 

l^iit t)I det g = + e{a) + ( J^^ + JVown) ^(-«) (72) 

In the unwarped, infinite flux case where a consise geometric result is obtained, one can in- 
tegrate out the axio-dilaton to obtain an effective density only in terms of the complex moduli. 
However, in our case this type of integration proves intractable. As a result we will when compar- 
ing with simulations have to fix a value of the axio-dilaton and compare the un-integrated form 
of our density. 



4 Numerical Vacuum Statistics 

To perform a numerical study of the distribution of vacua in moduli space near the conifold point, 
we will randomly choose appropriate fluxes F = {Fq, Fi, F2, F^) and H = {Hq, Hi, H2, H^) and 
then solve the conditions D^W = and D^^W = for the moduli space coordinate ^. Here 
W — Njlli is the superpotential, and N is an 8- vector whose first four components are those of F 
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and last four are those of H. We work in a basis such that the vector of 4-fold periods 11 = (E, tE) , 
where E is the vector of periods on the 3-fold. Near the conifold point, the vector of 3-fold periods 
takes the form: 

CXD 

S = J]a„r + &eiogHO, (73) 

n=0 

where the o„ and h are constant vectors associated with the expansion of the periods. Note that 
in the case of a single complex modulus, the vector h — (0, 0, 0, 6°), since only Eq has non-trivial 
logarithmic behavior near the conifold. Also, the local coordinate around the conifold point is 
proportional to E3, which implies that the vector ao = (0, Cq, Cq, Cq). 

4.1 Unwarped Analysis 

The unwarped Kahler potential is 

e-^ = -iS • g • E = -I ((a„ • Q ■ a^)Cr + {b ■ Q ■ a^)rC log(^0 + {a^-Q- h)Ci^og{-ii)) , 

(74) 

where the term proportional to b ■ Q ■ b has been dropped since given b and f] it vanishes. 
For SUSY vacua in the unwarped case 

D^W = N ■ (d^U + UK^) = 0. (75) 

Keeping logarithmic and constant terms gives 

(F - tH) • (ai + 6 {log{-zO + 1)) -{F- tH) ■ ao !° ' ^ ' = 0. (76) 

ao • V • Oo 

where the fact that clq - Q - h — Q has been used to simplify the expression. This is an equation of 
the form 

^ + Blog(-iO = 0, (77) 

with 

A = -{F- tH) ■ b{ao ■Q-ao) + {F- tH) ■ ai(ao ■ Q ■ a^) - {F - tH) ■ ao(ai • Q • Oo) (78) 
c 

B = {F-TH)-b{ao-Q-ao). (79) 
The leading-order constraints arising from requiring D^W — are 

r=^ = ^. (80) 



This imphes that 



H • oo 
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Before considering the effects of warp corrections, it's worth determining how close to the 
conifold vacua may be found in the unwarped scenario. The D^W = constraint imphes that 
|c^| is exponentially suppressed by the ratio of so if |^| is even just a couple of orders 

of magnitude greater than \B\, we should expect to see vacua on the order of 10^^^^ units away 
from the conifold point — indeed, this has been observed in previous studies. In order for |^| to 
differ appreciably from \B\ the quantity \{F — tH) ■ b\ should be relatively small compared to 
\{F — tH) ■ ao| or \{F — tH) - ail. Using the form of the vector b above, this indicates that the 
fluxes through the collapsing cycle, F3 and H3 should be small relative to some of the other fluxes. 



4.2 Warped Analysis 

Introducing warping leads to the corrections fH8|) and f H9|) to the Kahler potential and its deriva- 
tive. The modification to the near-conifold SUSY vacuum condition is then 



N ■ n. (82) 



Now, assuming that is small (i.e. the volume of the 3-fold is large) these new terms will matter 
only close to ,^ = 0. The SUSY condition thus leads to 

Fl/3 

A + B\ogi-tO+C^ = 0, (83) 

with A and B as before and 

C = 3C^{F -TH)-ao. (84) 

From this, we can see the rough influence of warping on the distribution of vacua. In the 
unwarped case, we expect to find vacua 10^^*^° or so away from the conifold with fluxes yielding 
|v4| ~ 100|S| (which with fluxes constrained to lie in (0, 100 is about the maximum order of 
magnitude difference that we expect.) If however, Cy^ ~ 10^^°, then for |^| ~ lO^^'^'', the warp 
term contribution is on the order of 10^°, swamping the logarithmic contribution and requiring 
fluxes ~ 10^° which lies beyond the range we consider. 

In the region of strong warping where the logarithmic term is dominated by the warping term, 
the distance of a vacuum from the conifold point is thus set by Given that A is at 

maximum of roughly 100 or so, the constant C^,, and thus, the overall volume of the Calabi-Yau, 
determines how near the conifold vacua lie. This can dramatically truncate the range — since the 
assumption of large but finite volume is well satisfied by volumes of order 10^°, but in those cases, 
vacua will not show up much closer than 10~^°. We can get vacua at around 10~^^° by taking a 
volume of order 10^*^, but in the absence of warping, vacua as far in as 10~^°° are expected. Thus, 
warping pushes vacua away from the conifold point. 
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4.3 Monte-Carlo vacua 



For the numerical analysis, we use the Calabi-Yau manifold labeled model 3 in the appendix of 
[1]. This family of Calabi-Yau can be expressed as a locus of octic polynomials in WP^'^'^'^'^. 
The corresponding orientifold arises from a certain limit of F-theory compactified on a Calabi- 
Yau fourfold hypersurface in WP^^'*^'^'^'^'^, following the methods of [lOj, and briefly described 
in [11]. For our purposes, we use the fact that the fourfold has Euler characteristic x = 23328, 
which implies that Lmax = x/24 = 972 for the tadpole condition for flux compactification on the 
corresponding orientifolded 3-fold. 

Since the warped form of the near conifold equation is not as simple to solve as in the unwarped 
case, a slightly more involved approach is necessary. We begin by defining two real variables p 
and 6 such that 

-i^ = p3e^^ (85) 

We take p > and < 6* < 27r. In terms of these variables, eqn (15^ and its complex conjugate 
expression take the form 

A + 'iB\n{p)+iBe + -^e-'^ = (86) 

~A + 3S ln(p) - m + -e'^ = (87) 

P 

Multiplying the first equation by B and the second one by i3, and then adding and subtracting 
the two, we find two purely real or imaginary equations. Letting A = ae*", B = 6e*'^, and C = ce^'^, 
we have 

asm{a-(3) + be + -sm{-f-l3-9) = (88) 
P 

acos(a-/3) + 361og(p) + -cos(7-/3-e) = (89) 

P 

We now solve for p in terms of 6 and then numerically solve the final equation for 6. It seems 
natural to solve equation ( 188|) for p since it is a linear equation. However, this approach fails in 
the limit C^, — )■ since then c — t- too. Instead we solve for p in equation fl89p . One can rearrange 
the equation as 

r 

pe^ \og{pe^) = -— cos(7 - /3 - 0) (90) 
3b 

Here we have defined the constant F = !iE2B^^z£L^ This is of the form xlog(a:;) = y which has 
the solution x = y/W{y) where W{y) is the Lambert ly-function. We therefore find 

Pie) = (91) 

3bW{-^cosi-f- (3-9)) 
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Vacuum Count 




Figure 2: A comparison between numerical and analytical distributions. Red circles mark the 
numerical data while the blue curve is the integrated analytical distribution. Distance from the 
conifold 1^1 is plotted on a log scale on the horizontal axis, while the vacuum count is plotted on 
the vertical axis. 

Consider equation (155]) . which now only depends on 6. Under the assumption that there is only 
one near conifold vacuum for each set of fluxes, the left hand side must either start out positive, 
and go negative or vice versa. To find the zero-crossing, we divide the region [0, 2tx\ into two 
equally pieces and then determine in which region (if any) equation (188 p changes sign. If such a 
region is found, we apply the same method to that region, splitting it into two smaller intervals, 
continuing in this way until we reach a predetermined level of accuracy. There are two relevant 
comments. First, in equation (19T]) it is not clear that the value of p is real, or even positive. We 
must therefore exclude the regions where p is either negative or complex. Fortunately, if p is real, 
it is never negative since W{x) must have the same sign as x. A necessary and sufficient condition 
for p to be real is that the argument of the Lambert W function is greater than or equal to — 1/e. 
This means that the relevant region to begin with may not be the entire interval [0, 27r]. Second, 
it turns out that the Lambert W function has two real branches for arguments between —1/e and 
0. Thus, both of these branches must be considered. 

To better compare the numerical and analytical and numerical distributions, we fix r and then 
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select a random sets of fluxes F and H consistent with our choice of r and satisfying the tadpole 
condition, F ■ Q ■ H < L^ax- For the particular model we consider, L^ax = 972, and we display 
a run using t = 2i, and Cy^ = 10~^^ in figure [2l The figure shows a numerical run compared to 
the analytical distribution. We plot the vacuum count and integrated analytical distribution as 
measured around the conifold point using a log scale for the distance from the conifold. As is 
evident from the figure, the count receives two major contributions: the one farther away from 
the conifold point is the usual contribution that is present without warping. However, we also 
see a major contribution much closer to the conifold at a distance roughly on the order of C^. 
This contribution is due to the strong warping effects and is matched by the cumulative analj^ical 
results. 

5 Discussion 

We've analyzed the distribution of flux vacua in the vicinity of the conifold point, including 
the effects of warping, and confirmed our results by a direct numerical Monte Carlo search. In 
comparison with the well known results, that don't include warping, we find a significant dilution 
of vacua in close proximity to the conifold, with the proximity scale set by the volume of the 
Calabi-Yau compactification. 

One complication in the analytical approach, relative to the unwarped case, is the need to 
bound the fluxes - a physically sensible requirement but one that can be avoided in the unwrapped 
analysis, yielding the geometrical result of [3lll]. It would be interesting to see whether the warped 
distribution of vacua can once again be related to intrinsic properties of the moduli space through 
a more complete geometrical treatment, likely requiringing careful consideration of the generalized 
complex geometry of conformal Calabi-Yau spaces [HIIZ]- 
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A Appendix 



A.l Covariant Derivatives 

We start with the standard definition of the Kahler potential 



A rescahng of the holomorphic 4-form Q4 — > €^^^^0,4 imphes that K ^ K — f(z) — f{z). The 
covariant derivative is defined so that it is covariant under such rescahngs: 



implying that DoVLa = {da + Ka) f^4. Note also that the holomorphic covariant derivative annihi- 
lates antiholomorphic objects, i.e. D0O4 — da^4 — 0. We see then that 

Dae-"" = {da + Ka) e"^ = 

Notice that covariance dictates that 

Dae"" ^{da-Ka) 6^ = 

In addition to this Kahler scaling structure, the complex structure moduli space is a manifold 
with a natural Kahler metric K^^ = dad^K. We can thus define a metric compatible connection 
via 

where K^^^ — dcK^. Note that the connection components with mixed holomorphic and anti- 
holomorphic indices vanish. By suitably extending the covariant derivative we can ensure 
that it transforms covariantly under both Kahler rescalings and coordinate transformations on the 
complex moduli space. In particular, since the connection is metric compatible, DaK^c — 0. 
The superpotential 

W ^ j G4/\n4 

scales as Q4, and thus 

DaW^daW + KaW 

A super symmetric vacuum satisfies the conditions DaW — 0. The components of the fermion 
mass matrix are 

daD.W = dadbW + KabW + K,Wa 
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Notice that at a generic point in the moduh space the quantity 

DaD.W = daD.W + KaD.W - Tl,DM = daD.W + KaD.W - K,,jK^'dJv 

does not equate to the fermion mass matrix. However, the extra terms drop out at supersymmetric 
vacua. 

Now let ^4 — )■ ^4 = e^/^f24 and similarly for the (0,4)-form. Notice that the scaling properties 
of ^4 imply that 

DaQ^ = (^da + ^Ka^ ^4 = (^0^ + ^K^^ (e^'/^Q^^ = e''/' (9„ + K^) = e'^'^DSt^ 

The rescaled (4,0) and (0,4) forms are convenient since they remove factors of from various 
expressions. In particular we have 

j ^4 A ^4 = 1 

We can also go to an orthonormal frame by introducing vielbeins Sj^j^ = e\e^K^i. The covariant 
derivative must be extended so as to keep the vielbeins covariantly constant: 

implying that 

Given these definitions, we can now go to rescaled expressions in the orthonormal frame: 

DaDbW = OaDbW + KaDbW - ujab^DcW 

once again, the expression above agrees with the fermion mass matrix components evaluated at a 
vacuum. 

A. 2 Computations 

In this section, we provide derivations for the equations (13T]) -( 138|) . To do so, recall from equation 
( 130|) that the flux four form written in the basis B = Da^a, DqDjQ^, ^4, Da^a, DqDjQ^} is 

Ga = Xn^ - y^DA^i + Z^DoDjn^ + c.c. (92) 

It's useful to note that DaDb^a is a (2,2)-form. In fact, given the nature of our Calabi- 
Yau 4-fold, essentially factorizing into a 3-fold and a torus, this (2,2)-form can be decomposed 
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as (2, 2) = (1, 0) A (1, 2) ® (0, 1) A (2, 1). To see this, note that 0^0,4 could be a mixture of a 
(3, 1) © (4, 0), but the (4,0) component vanishes: 

/ Da^4 a n4 = Da ( / A ^4 ) - / O4 A Da^4 = 
Jm \Jm J Jm 

where the two terms vanish given the properties of the covariant derivative defined above. Similarly 
DaDb^4 could in principle have (2, 2) ® (3, 1) © (4, 0) structure. However, 

/ DADB^4AUi = Da( Db^^a'^A- Db^4ADa^4 = 
Jm \Jm J Jm 

implying that there is no (4, 0) component. Furthermore 

/ DaDb^a a D^n4 = Da( Db^4 a Dc^4 ] - Db^4 a DaD^'^4 
Jm \Jm J Jm 

the first term on the left-hand-side is equal to DaSbc — 0- The second term becomes 

^AC / Db^4 = 

Jm 

which shows that there is no (3, 1) component in DaDb^4- We now turn to the identities of 
interest. 

m W = X 

By the definition of the superpotential, we have 

W = /" G'4AQ4 
Jm 

= / {Xfl4 - Y^DAfl4 + Z^DoDin4 + c.c.) A O4 = ^ (93) 
Jm 

In the last step we used the orthonormality of the basis. 
. DaW = Ya 

Once again, we will use the orthonormality of the basis. In particular we have (since G4 is 
independent of the moduli) 



DaW = / G4ADAn4 
Jm 

= / {Xn4 - F^L'Afi4 + Z^DoDiQ4 + c.c.) A DAn4 
Jm 

= -Y^ [ DbT14 a Da^4 = +Ya (94) 
Jm 

In the last step we again used (J^ Db^4 A Da^4 — Sba)- 
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• DoDqW = 
We have 

DoDqW = [ GihDoDo^i (95) 
Now D^4^ is a (0, 1) A (3, 0)-form. In fact, we see that 

Dr^i = {Or + Kr) {a - t/3) = Kr {a - T/3) = KrVti 

where we have used — —1/ [t — t). Using the fact that the vielbein eg = l/-f^T, we have 
-Dgf^4 = Qi A Q3, however we know that -Dq^i = 0, so the identity holds. 

• DqDiW = Zi 

This identity follows from orthonormality: 

D^DiW = [ G4ADoDin= [ {Xn^-Y^DAn^ + Z^DoDin^ + cc.) ADoDin = Zj 
Jm Jm 

(96) 

. DjDjW = TukZ"^ 

We will again use the definition for G4 

DjDjW = [ (XQi - F^D^fi4 + Z^DoDiQi + c.c.) A DiDjQ^ (97) 
Jm 

As discussed above D^Ds^i is a (2,2)-form which breaks up as (1, 0) A (1, 2) © (0, 1) A (2, 1). 
Now, DiDjVt^ is precisely a (1, 0) A (1, 2)-form, so the covariant derivatives only act on the 
3-fold factor. The only pieces of the integral above that can yield a non-zero result must be 
of the form (0, 1) A (2, 1), which are thus proportional to DoDifl^. This leaves us with 



DjDjW = Z^ [ DoDk^^ a DjDjQ^ 
Jm 

The Dk and Dq derivatives commute, so we have 

DjDjW = Z^Dk ( [ DoQi A DiDjQA -Z^ [ D^Qi A DkDiDjQ^ 
\Jm ) Jm 

The first term on the right-hand-side vanishes due to orthonormality since DqQ^ is a (3,1)- 
form while DjDjVli is a (1, 0) A (1, 2)-form. Factorizing = fii A ils, we see that ^0^4 = 
A O3. The integral over the torus will simply yield a factor of — i, leaving us with 

DiDjW = iZ^ f 1^3 A DkDjDjQ^ 



cy 
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However, pulling out all scaling factors and vielbeins, we see that the resulting derivatives 
can all be converted to partials. This allows us to rearrange the ordering and gives 

DiDjW = iZ"^ [ nsA DjOjOk^s = TukZ"^ 
Jcy 

• DaDbW = 5abX 
First consider 

The first and last term vanish since W is holomorphic in the moduli. Then, since W = X, 
by reintroducing the scaling factor e^^^ and the vielbeins, we have 

DaDbW = 6abX (98) 

• DqDiW = 

We can easily see this by noting that the outer derivative is a regular partial derivative and 
that this commutes with the inner derivative. Then, since W is holomorphic in r, Bq sends 
the expression to zero. 
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